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Abstract: We develop a general method of calculating Stimulated 
Brillouin Scattering (SBS) gain coefficient in axially periodic waveguides. 
Applying this method to a silicon periodic waveguide suspended in air, 
we demonstrate that SBS nonlinearity can be dramatically enhanced 
at the brillouin zone boundary where the decreased group velocity of 
light magnifies photon-phonon interaction. In addition, we show that the 
symmetry plane perpendicular to the propagation axis plays an important 
role in both forward and backward SBS processes. In forward SBS, only 
elastic modes which are even about this plane are excitable. In back- 
ward SBS, the SBS gain coefficients of elastic modes approach to either 
infinity or constants, depending on their symmetry about this plane alq = Q. 
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1. Introduction 

Stimulated Brillouin Scattering (SBS) is a third order nonlinear process in which two opti- 
cal modes are coupled through an elastic mode ||T[|2]. In a waveguide system, the interference 
of pump and Stokes waves generates dynamic optical forces at the beat frequency. The opti- 
cal force, while on resonance with an elastic mode at the phase-matching wavevector, excites 
the mechanical vibration in the waveguide, which in turn scatters light between the pump and 
Stokes waves. Since its discovery, SBS has been extensively studied with a variety of applica- 
tions in efficient phonon generation [T'4], optical frequency conversion fS'-'Tl, slow light ||8]-[TTj 
and signal processing techniques | ,12,,13J . 

Previous experimental and theoretical studies of SBS have heavily focused on axially in- 
variant systems, such as optical fibers, phononic crystal fibers, and rectangular waveguides 
l@HZl[llhlZJ' desirable to study the SBS process in axially periodic waveguide for two 
reasons. First, the unit cell in periodic structures provides more degree of freedoms to tailor the 
distributions of optical force and elastic deformation fVf\. More importantly, the slow group ve- 
locities of light, which are readily achievable in periodic waveguides, can dramatically enhance 
SBS nonlinearities [18]. In axially invariant waveguides, all the vectorial fields can be decom- 
posed into real/imaginary-valued longitudinal and transverse components. Such decomposition 
are useful in characterizing various SBS processes |17|. In axially periodic waveguides, all the 
components of electric field and elastic deformations become complex-valued. A careful anal- 
ysis of the symmetry in optical forces and elastic modes is required to elucidate photon-phonon 
coupling and the resultant selection rules in various SBS processes. 

In this article, we propose a general method of calculating SBS gains in periodic waveguides. 
Armed with this general formalism, we study the forward SBS (FSBS) and backward SBS 
(BSBS) processes of a suspended silicon periodic waveguide. We demonstrate that the strength 
of both FSBS and BSBS can be dramatically enhanced in the slow light regime. In addition, we 
show that the symmetry plane perpendicular to the propagation axis plays an important role in 
characterizing different SBS processes. This plane separates elastic modes at ^ = into even 
and odd modes. For FSBS, only even modes are excitable. For BSBS, the gain coefficients 
approaches to either infinity or constants, depending on the mode symmetries s.lq = Q. 



2. Calculating the SBS gain of periodic waveguides 



To start with, we develop a general method of calculating the SBS gain of periodic waveguide 
systems. The optical (elastic) eigen-modes of a waveguide with axial periodicity a is char- 
acterized by Bloch wavevector k (q). In a typical SBS process, the pump wave Epe^"''p' at 
wavevector kp and Stokes wave E^e^""^' at wavevector k^ generate optical forces varying in 
space at wavevector q = kp — k^ and oscillating in time at the beat frequency Q. = cOp — Ov 
This optical force can excite mechanical vibrations which enables the parametric conversion 
between pump and Stokes waves. Assuming the propagation direction is parallel to the x axis, 
this process can be described by the following relation (T) 



dx 



gPpPs-asPs (1) 



Here, Pp and P^ are the guided power of the pump and Stokes waves, and g is the SBS gain. 
Through particle flux conservation, SBS gain is given by the follow formula |17) : 

s{^) = J^ReU'^\ (2) 



laaPpPs \' dt / 

where f is the optical force generated by pump and Stokes waves, and u is the elastic response 
of the waveguide induced by f. The optical power of a periodic waveguide is given hy P = 
Vg(E,eE)/2fl, where v^, is the optical group velocity. Therefore, 

' y,pV,,(Ep,eEp)(E„eEs) 
The elastic response u can be decomposed into elastic eigen-modes at q, u = b,„u,n with b,„ 



given by Eq. ( A.ll i. The total SBS gain is the sum of SBS gains of individual elastic modes: 



The SBS gain of the m elastic mode has a Lorentian shape and a peak value G„,: 

^m^mVgpVgs (u,„pUm) (Ep, gEp) (Eg, gEj) 

In Eq. Q, the overlap integrals of different optical forces are linearly summed. Here, we con- 
sider three kinds of optical forces: electrostriction body force, electrostriction pressure, and 
radiation pressure | [T7][T9) . Electrostriction body force is integrated over the volume of the unit 
cell, while electrostriction pressure and radiation pressure are integrated over the boundary of 
the unit cell. Equation (jSj) shows that the SBS gain coefficient is dependent on the frequency 
ratio, the mechanical loss factor, the optical group velocities of pump and Stokes waves, and 
the overlap integral between optical forces and elastic eigen-modes. With periodic waveguides 
with unit cell structures, not only can the overlap integral be tailored, we can also take advan- 
tage of the decreased optical group velocity near the brillouin zone boundary to gain further 
enhancement of SBS nonlinearities. 



3. Optical and elastic modes of a silicon periodic waveguide 



Using the formalism developed above, we proceed to study the SBS process of a a suspended 
silicon waveguide with periodic cylindrical holes (Fig. [TJa) insert). The axial periodicity is a, 
the cross-section in yz plane is a by QAa, and the radius of the cylindrical air hole is 0.25a. 
For silicon, we use refractive index n ~ 3.5, Young's modulus E — 170 x 10^ Pa, Poisson's 
ratio V = 0.28, and density p ^ 2329kg/m2. In addition, we assume that the [100], [010], 
and [001] symmetry direction of this crystalline silicon coincide with the x, y, and z axis re- 
spectively. Under this orientation, the photo-elastic tensor piji^i in the contracted notation is 



[/7ii,/7i2,p44] = [—0.09,0.017,-0.051] |20|. The waveguide has three symmetry planes x = 0, 
J = 0, and z = 0. We assume that the crystalline structure of the material is also symmetric 
about these three planes so that the anisotropy in optical, elastic, and photo-elastic constants 
doesn't break these mirror symmetries. This condition is clearly satisfied for silicon in the cur- 
rent orientation. 

First, we analyze the optical modes of the waveguide. The optical modes are categorized 
into yeven/yodd (zeven/zodd) according to their symmetries about plane y = (z = 0). The 
fundamental mode is yodd and zeven with Ey as the dominant component of electric field (Fig. 
[TJa)). We fix the pump wavelength at 1.55/im. So a different operating point in the dispersion 
relation indicates a different value of a. The optical mode E doesn't have symmetry about 
plane x = since nonzero k breaks this mirror symmetry. Actually, the mirror reflection of 
eigen-mode E at ^ corresponds the eigen-mode at —k, which is also the complex conjugate of 
E. Under a properly chosen phase, the mirror reflection of E is exactly E*. We can write this 
relation as: 

Ei{-x,y,z)=E*{x,y,z)si (6) 

where S;^ = — 1 , and Sy, = 1 . The subscript in s does not introduce summation when encountered 
with repeated indices. 

Next, we analyze the elastic modes of the waveguide. Again, the elastic modes are catego- 
rized into yeven/yodd (zeven/zodd) based on their symmetries about y = (z = 0). In intra- 
modal coupling where both pump and Stokes waves reside in the same optical eigen-mode, the 
optical force is always symmetric with respect to planes y — and z = Q- Therefore, we only 
need to consider E-modes which are both yzeven and zeven (Fig.[ljb)). At g' = 0, the symmetry 
about plane x = is recovered, separating E-modes into xeven and xodd modes (Fig.[TJc)). In 
addition, the elastic eigen-equation at g = is invariant under conjugation operation, resulting 
in real-valued u. At nonzero q, the elastic modes are neither symmetric nor anti-symmetric 
about plane x = 0, and u are complex-valued. Similar to the optical modes, we can choose a 
proper phase of u so that the mirror reflection of u is exactly u: 

Ui{-x,y,z) = u'^{x,y,z)si (7) 

4. Forward SBS 

In FSBS, pump and Stokes waves approximately correspond to the same optical mode Ee^'™, 
and excite standing-wave elastic modes atq — 0. Under such conditions, electrostriction tensor 



in Eq. (B.l i and MST in Eq. (B.5 1 are given by: 



CT.j - -^eon*pijk,Re{EkEl) (8) 
T,j = Uoei2ReiEiE*)-5ij\E\^) (9) 



Both electrostriction tensor and MST are real-values, resulting in real-valued optical forces. 
We select an operating point at O) = 0.265{2nc/a) and k — 0.15{n/a) with a — 420nm, and 
compute the distributions of electrostriction body force, electrostriction pressure, and radiation 
pressure (Fig.|2];a)). The dominant component of electrostriction body force is ff^, because the 
dominant component of electric field is Ey, and pi i is about five times larger than pi2- Radiation 
pressure point outwards, which is about five times larger than electrostriction pressure. 

One important feature about the optical force distribution is that all optical forces are sym- 
metric about plane x = Q although the optical eigen-mode E doesn't have this symmetry. This 
can be formerly proven by examining the symmetry of electrostriction tensor and MST. Be- 
cause the crystal structure of the waveguide material is symmetric about x = 0, the photo-elastic 
tensor is zero when there is odd number of x in the subscript: piji^i — PijkiSiSjSkSi- Using this 
property and = we have 

Oij{-x,y,z) = -^eon'^PijkiRe{Ek{-x,y,z)El{~x,y,z)) 

= -^eon'^PijkiRe{El {x, y, z)Ei {x,y, z) jSkSiSkS, 

= <^ij{x,y,z)siSj (10) 

Similarly, using = we have 

Tij{-x,y,z) = ^eoe{2Re{Eii-x,y,z)E*{-x,y,z)) - 5ij\E{-x,y,z)\^) 

- ^eoe{2Re{E,{x,y,z)E*{x,y,z)) - 5ij\E{x,y,z)\^)siSj 

= Tij(x,y,z)siSj (11) 



Combining Eq. ( 10 1 and (111 with the fact that optical force is given by the divergence of the 
corresponding tensor, we conclude that both optical forces in FSBS are symmetric about plane 
x^O. 

The symmetry of optical forces, together with symmetry property of elastic modes at <7 = 0, 
indicates that only xeven modes are excitable. We calculate the FSBS gain assuming a mechan- 
ical quality factor Q = 1000 for all the elastic modes. As expected, only xeven modes E2 and 
E4 have nonzero FSBS gains (Fig. |2|b)). Mode E2 has large displacement in y direction and 
small displacement in z direction. Such modal profile agrees well with electrostriction body 
force and the radiation pressure on the outer lateral surfaces, generating large FSBS gains from 
electrostriction (0.71 x lO^m'^W"') and radiation pressure (0.44 x 10^m"'W"'). Furthermore, 
these two effect add up constructively, resulting a total FSBS gain as large as 2.27 x lO'^m"' 
Mode E4 has a small FSBS gain of 0.17 x lO^m 'W"' because (1) the nodal planes of Uy reduce 
the overlap integral of electrostriction body force and (2) radiation pressures on the outer and 
inner lateral surfaces are canceled out to a large extent. 

Next, we study how the FSBS gain varies as the operating point moves from brillouin zone 
interior to boundary (Fig.[3]l. For mode E2, electrostriction force and radiation pressure always 
add up constructively, creating an even larger FSBS gain. For mode E4, the radiation-pressure- 
only gain coefficient vanishes around k = Q.51{n/a), because of the cancellation of radiation 
pressures on different surfaces. For both mode E2 and E4, when the operating point approaches 
brillouin zone boundary, the overlap integrals approach to constants while the optical group 
velocity vanishes as 0{Ak) (Ni ^\k — n/a\). As a result, the FSBS gains approach to infinity 
as 0{\/Ak^). 



5. Backward SBS 



In BSBS, pump and Stokes waves travel in the opposite directions, exciting elastic modes at 



q = 2k. Under such conditions, electrostriction tensor in Eq. (B.l i and MST in Eq. (B.5 i are 
given by: 

o-/,/ = -^eon'^PijkiEkEi (12) 

T,j = ^eoe{2E,Ej-5ijE-E) (13) 

Both the electrostriction tensor and MST are complex-valued, resulting in complex-valued 
optical forces. We select an operating point at o = 0.265{2nc/a) and k — Q.15{n/a) with 
a = 420nm, and compute the real and imaginary parts of electrostriction body force, elec- 
trostriction pressure, and radiation pressure (Fig. |4|a)). Similar to the case of FSBS, we can 
show that the electrostriction tensor and MST in BSBS have the following properties: 

(yij{-x,y,z) = a*j{x,y,z)siSj (14) 
Tij{-x,y,z) = T*j{x,y,z)siSj (15) 

Taking the real and imaginary parts of the expressions above, it is straightforward to show 
that the real (imaginary) part of optical forces is symmetric (anti-symmetric) about plane x = 
0. Under such optical forces, all the elastic modes at q ~ 2k are excitable. We calculate the 
BSBS gain assuming a mechanical quality factor Q = 1000 for all the elastic modes (Fig.|4|b)). 
Mode E2 has the largest BSBS gain (1.42 X lO'^m-'W-i), which comes from a constructive 
combination of electrostriction (0.29 x 10'*m"^W"') and radiation pressure (0.42 x 10'*m"' W"'). 

One interesting feature about the gain coefficients is that electrostriction and radiation pres- 
sure add up either purely constructively or destructively: 

(16) 

This can also be explained using symmetry argument. As mentioned above, the real (imaginary) 
parts of optical forces are symmetric (anti-symmetric) about plane x — 0. On the other hand, 
the real (imaginary) part of elastic displacement is even (odd) about plane x = 0. Although both 
optical forces and elastic modes are complex-valued, their overlap integrals are always real. The 
direct interference between different optical forces can be exploited to deliberately enhance or 
suppress the SBS nonlinearity of certain elastic modes. 

Next, we study how the BSBS gain varies as the operating point approaches the slow light 
regime (Fig.|5]). For mode E2, electrostriction force and radiation pressure always add up con- 
structively. In the slow light regime, the BSBS gain approaches to infinity as 0(1 /A^^) as the 
optical group velocity vanishes. In contrast, the BSBS gain of mode E3 approaches to a con- 
stant at the brillouin zone boundary. This comes from two properties associated with k = nja. 
At A; = %la, the optical modes at k and —k merge into one mode, and BSBS becomes equiva- 
lent to FSBS. So the optical forces in BSBS become symmetric about plane x = 0. For k close 
to nja, the optical force in BSBS can still be decomposed into symmetric and anti-symmetric 
components with the anti-symmetric component on the order of OiJSk). On the other hand, at 
k = n / a, q = 2k — 2n / a, which is equivalent to ^ = 0. So the elastic modes in BSBS recover 
their symmetries about x = Q. For k close to n/a, elastic modes can be decomposed into odd 
and even components with respect to plane x = 0. The even component in xodd modes such as 
mode E3 is on the order of (9(AA:). Therefore, in BSBS, the overlap integral between optical 
forces and mode E3 vanishes as (9(AA:). As a result, the BSBS gain approaches to a constant 



rather than infinity in the slow fight regime. 
6. Concluding remarks 

In tfiis article, we analyze tfie forward and backward SBS processes of a periodic waveguide 
suspended in air. Tfie suspended structure provides tigfit lateral confinement of ligfit and nearly 
perfect lateral confinement of sound p7) . Tfie periodic structure slows down tfie optical group 
velocity. Tfie combination of tfiese two effects creates a giant enfiancement of SBS nonlin- 
earity over conventional nonlinear fibers | [T4[ . In addition, we cfiaracterize elastic modes and 
tfie resultant SBS gain coefficients based tfieir symmetries witfi respect to tfie symmetry plane 
perpendicular to tfie propagation axis. Our analysis doesn't rely on tfie specific waveguide ge- 
ometry or tfie crystalline structure of silicon. Tfie conclusion about tfie relation between mode 
symmetry and SBS gain is valid as long as (1) tfie waveguide fias a symmetry plane perpendic- 
ular to tfie propagation axis and (2) tfie crystalline structure of tfie underlying material doesn't 
break tfiis symmetry. Our analysis can be readily applied to simpler structures sucfi as axially 
invariant waveguides and more complicated structures sucfi as waveguide systems witfi pfio- 
tonic/pfiononic lateral confinement. 



A. Eigen-mode Decomposition of Elastic Response 

A.l. Orthogonality of elastic eigen-modes 



When loss and external forces are ignored, the equation about displacement u is 1 21 1: 

where p is the mass density, and c,y<./ is the elasticity tensor. Cjji^i is symmetric about the first 
two and last two indices: Cjji^i — Cji^i, ciji^i — Cijik. It is also symmetric when the the first two 
indices and the last two indices are interchanged: c/tUj ~ ciji^i \TP\. 

For a finite structure or an infinite structure with periodicity, the elastic eigen-modes are 
discrete. Considering an elastic eigen-mode with frequency i2: u = ue^'^', we have: 



-Cijki^ = -Q.yui (A.2) 
ax ! oxl- 



This is an eigen-equation with operator F and kernel p : 

Fn ^ -il^pu (A.3) 

We define the inner product between two vector fields as the overlap integral over the finite 
structure or the unit cell of an infinite structure with periodicity: 

(A,B) = jA*BidV (A A) 

By integrating by parts and using the properties of Cjji^i, we can show that operator F is Hermi- 
tian: 

d dBi 

' dxj '^^^ dxji 

dA* dBi 
Cijki^ — -^aV 
axj axji 



(A,FB) = lA*—Ciju^dV 



dBi dA; 

Cijkl^^ — 
axj axji 

d dAI 
Bi^Cijki^dV 
axj axji 

= (FA,B) (A.5) 
By the property of Hermitian operators, elastic eigen-modes are orthogonal under kernel p : 

(u,„,pu„) = 5„,„(um,pu,„) (A.6) 
A.2. Elastic response with loss and external forces 

Now we apply external force f to the structure. The elastic equation becomes: 



Assuming / is time-harmonic with frequency we have 

^Cijki jr^+fi- -p^\ (A.8) 
axj axii 

By decomposing u into the eigen-modes u = Lm^;nU„„ we get; 

^/7,„p(£l2 „n2)u,„ = f (A.9) 

m 

Taking the inner product of u,,, with both sides of the equation above and applying the orthog- 
onaUty condition, we have: 

(u,„,pum) n^-n-^ 

We now consider a more general case, where the mechanical loss is present. The commonly 
encountered mechanical loss mechanism includes air damping, thermoelastic dissipation, and 
clamping losses [22 J . The first order effect of loss can be captured by adding an imaginary part 
to fi„,: n,„ — ir,n/2. The mechanical quality factor is defined as Q,„ = i2,„/r,„. Assuming Q,„ 
is well above 1 , we have, 

u _ (um,f) 1 r„,/2 



B. Calculation of Optical Forces 

B.l. Electrostriction forces 

Electrostriction force is derived from electrostriction tensor. The instantaneous electrostriction 
tensor is given by: 

1 4 

(7ij = --eon pijkiEkEi (B.l) 

where n is the refractive index, and pijiii is the photoelastic tensor |23|. When both pump and 
Stokes waves are present, the total electric field is given by [Epe^""''' +Ese^"'''')/2 + c.c. 
Inserting this expression to Eq. (B.l i, and taking out the components with frequency £2, we get 
the time -harmonic electrostriction tensor aije^'^': 

Gij - -^-em^pijk,{EpkE:,+E„,E*^) (B.2) 

Electrostriction force is the divergence of electrostriction tensor. In a system consisting of 
domains of homogeneous materials, electrostriction forces can exist within each material (elec- 
trostriction body force) and on the interfaces (electrostriction pressure). Electrostriction body 
force is given by: 

ff' = -djOij (B.3) 

Electrostriction pressure on the interface between material 1 and 2 is given by (normal vector n 
points from 1 to 2): 

FI"^ = {auj - (T2,;)n./ (B.4) 



Equation (B.4i shows that electrostriction pressure can have tangent components on material 
boundaries. 



B.2. Radiation pressure 

Radiation pressure is derived from Maxwell Stress Tensor (MST). For a dielectric system (ji ~ 
1) without free charges (p = 0,7 = 0), radiation pressure is localized where the gradient of e 
is nonzero |24 25). For a system consisting of homogeneous materials, radiation pressure only 



exists on the interfaces. The electric part of instantaneous MST is; 

Tij = eoeiEiEj-^-dijE^) (B.5) 
When both pump and Stokes waves are present, the total electric field is given by (Epg^""''' + 



E^e ""'') /2 + c.c. Inserting this expression into Eq. (B.5 > and filtering out the component with 
frequency H, we arrive at the time-harmonic MST Tije^^': 

T,i = ^-eoe{E,„E:j+E„jE:, - d,jE„ ■ E,*) (B.6) 

Radiation pressure on the interface between material 1 and 2 is given by the discontinuity of 
MST: 

Ff"' ^{T2ij-Tuj)nj (B.7) 

By decomposing the electric field into its normal and tangent components with respect to the 
dielectric interfaces E = E„n + Eti, and using the boundary condition e\Ein — £2E2n — Dn and 
Eit = E2t = Et, we can get the time-harmonic radiation pressure F^'^e^'^'': 

- -\eQE,,,E:,{e2-£i)n+U^'Dp„D*„{e^'-£-')n (B.8) 



Equation ( |B.8[ ) shows that radiation pressure is always normal to the interface. 




k{7r/a) ' qiir/a) 
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Fig. 1. The optical and elastic modes of a silicon periodic waveguide, (a) The dispersion 
relation of the fundamental optical mode which is yodd and zeven. The area shaded in gray 
represents light cone of air. (b) The dispersion relations of elastic modes which are yeven 
and zeven. Xeven and xodd modes at g = are colored in red and blue respectively, (c) 
Elastic modal profiles at ^ = 0. The deformation is proportional to the displacement u. The 
colored surface represents displacement component iiy. Red, white and blue correspond to 
positive, zero, and negative values respectively. Mode El experiences a parallel shift along 
X direction. 




Fig. 2. Optical force distributions and gain coefficients of forward SBS. Tire operating 
point is CO = 0.265(2;rc/a) and k = 0.75(;r/a) with a = 420nm. Elastic modes at g = are 
excited, (a) Distributions of optical forces. For electrostriction body force, the two subplots 
show the body force density on planes z = and y = ±0.4a. For electrostriction pressure 
and radiation pressure, the two subplots show the pressure on the lateral surfaces and the 
top surface. Electrostriction pressure is multiplied by 5 so that it can be plotted on the same 
scale as radiation pressure. All the optical forces are symmetric about to plane x = 0. (b) 
The FSBS gains of individual elastic modes assuming Q = 1000. Blue, green, and red bars 
represent the FSBS gains under three conditions: electrostriction-only, radiation-pressure- 
only, and the combined effects. 




Fig. 3. Forward SBS gains of mode E2 and mode E4 as k varies from Q.5{'Jt/a) to Tt/a, and 
the corresponding a varies from 354nm to 434nm. For botli mode E2 and E4, tlie FSBS 
gains approacli to infinity as 0{l/ M^) at thie brillouin zone boundary. 
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Fig. 4. Optical force distributions and gain coefficients of backward SBS. The oper- 
ating point is ta = 0.265(2;rc/a) and k = 0J5{n/a) with a = 420nm. Elastic modes at 
q= 1.5(;r/a) are excited, (a) Distributions of the real and imaginary parts of optical forces. 
For electrostriction body force, the two subplots show the body force density on planes 
z = and y = ±0.4a. For electrostriction pressure and radiation pressure, the two subplots 
show the pressure on the lateral surfaces and the top surface. Electrostriction pressure is 
multiplied by 5 so that it can be plotted in the same scale as radiation pressure. The real 
part of optical forces is symmetric about plane x = 0, while the imaginary part of opti- 
cal forces is anti-symmetric about plane x = 0. (b) The BSBS gains of individual elastic 
modes assuming Q = 1000. Blue, green, and red bars represent the BSBS gains under three 
conditions: electrostriction-only, radiation-pressure-only, and the combined effects. 




Fig. 5. Backward SBS gains of mode E2 and mode E3 as k varies from Q.5{lt/a) to n/a 
and the corresponding a varies from 354nm to 434nm. For mode E2, the BSBS gain ap- 
proaches to infinity as 0{i/hk^). For mode E3, the BSBS gain approaches to a constant at 
the brillouin zone boundary. 



